Information percolation and cutoff for the random‐cluster model

We consider the random‐cluster model (RCM) on with parameters p∈(0,1) and q ≥ 1. This is a generalization of the standard bond percolation (with edges open independently with probability p) which is biased by a factor q raised to the number of connected components. We study the well‐known Fortuin‐Kasteleyn (FK)‐dynamics on this model where the update at an edge depends on the global geometry of the system unlike the Glauber heat‐bath dynamics for spin systems, and prove that for all small enough p (depending on the dimension) and any q>1, the FK‐dynamics exhibits the cutoff phenomenon at with a window size , where λ_∞ is the large n limit of the spectral gap of the process. Our proof extends the information percolation framework of Lubetzky and Sly to the RCM and also relies on the arguments of Blanca and Sinclair who proved a sharp mixing time bound for the planar version. A key aspect of our proof is the analysis of the effect of a sequence of dependent (across time) Bernoulli percolations extracted from the graphical construction of the dynamics, on how information propagates.     

A proper orthogonal decomposition variational multiscale meshless interpolating element-free Galerkin method for incompressible magnetohydrodynamics flow

In the recent decade, the meshless methods have been handled for solving most of PDEs due to easiness of the meshless methods. One of the popular meshless methods is the element-free Galerkin (EFG) method that was first proposed for solving some problems in the solid mechanics. The test and trial functions of the EFG are based on the special basis. Recently, some modifications have been developed to improve the EFG method. One of these improvements is the variational multiscale EFG procedure. In the current article, the shape functions of interpolation moving least squares approximation have been applied to the variational multiscale EFG technique for solving the incompressible magnetohydrodynamics flow. In order to reduce the elapsed CPU time of simulation, we employ a reduced-order model based on the proper orthogonal decomposition technique. The current combination can be referred to as the reduced-order variational multiscale EFG technique. To illustrate the reduction in CPU time used as well as the efficiency of the proposed method, we applied it for the two-dimensional cases.     

Frequency-dependent acoustic energy focusing in hexagonal ceramic micro-scaffolds

Acoustic properties of an additive-manufactured SiC scaffold with hexagonal symmetry fabricated by the robocasting method are studied both numerically and experimentally. The numerical analysis is based on the finite element method (FEM) using Bloch boundary conditions. The calculations show both angular and frequency dispersion of the acoustic waves with wavelengths comparable to the spacing between the rods, i.e., on a millimeter scale, indicating interesting acoustic properties in the MHz range. The dispersion character leads to focusing of the energy propagation into the directions of the rods of the hexagonal structure. This is illustrated by modal-based calculations of the propagation of longitudinal and out-of-plane shear wave packets with a dominant wavelength. The experimental analysis consists of two steps, the measurement of the resonant spectrum and shear wave propagation character. The measured resonant spectrum is in good agreement with the one calculated using numerically obtained low-frequency properties of the structure, also showing the quality of the overall manufactured structure. The time-domain measurement shows significant changes in the energy propagation between low and high frequencies, as predicted by FEM calculations.     

Synthesis and characterization of VO–vanillin complex immobilized on MCM-41 and its facile catalytic application in the sulfoxidation reaction,and synthesis of 2,3-dihydroquinazolin-4(1H)-ones and disulfides in green media

In this work, a vanillin complex is immobilized onto MCM-41 and characterized by FT-IR, X-ray diffraction, scanning electron microscopy, energy dispersive spectroscopy, thermogravimetric analysis, and BET techniques. This supported Schiff base complex was found to be an efficient and recoverable catalyst for the chemoselective oxidation of sulfides into sulfoxides and thiols into their corresponding disulfides (using hydrogen peroxide as a green oxidant) and also a suitable catalyst for the preparation of 2,3-dihydroquinazolin-4(1H)-one derivatives in water at 90°C. Using this protocol, we show that a variety of disulfides, sulfoxides, and 2,3-dihydroquinazolin-4(1H)-one derivatives can be synthesized in green conditions. The catalyst can be recovered and recycled for further reactions without appreciable loss of catalytic performance.     

A new database and benchmark of the bond energies of noble-gas-containing molecules

We have developed a new database of structures and bond energies of 59 noble-gas-containing molecules. The structures were calculated by CCSD(T)/aug-cc-pVTZ methods and the bond energies were obtained using the CCSD(T)/complete basis set method. Many wavefunction-based and density functional theory methods have been benchmarked against the 59 accurate bond energies. Our results show that the MPW1B95, B2GP-PLYP, and DSD-BLYP functionals with the aug-cc-pVTZ basis set excel in predicting the bond energies of noble-gas molecules with mean unsigned errors (MUEs) of 2.0 to 2.1 kcal/mol. When combinations of Dunning's basis sets are used, the MPW1B95, B2GP-PLYP, DSD-BLYP, and BMK functionals give significantly lower MUEs of 1.6 to 1.9 kcal/mol. Doubly hybrid methods using B2GP-PLYP and DSD-BLYP functionals and MP2 calculation also provide satisfactory accuracy with MUEs of 1.4 to 1.5 kcal/mol. If the Ng bond energies and the total atomization energies of a group of 109 main-group molecules are considered at the same time, the MPW1B95/aug-cc-pVTZ single-level method (MUE = 2.7 kcal/mol) and the B2GP-PLYP and DSD-PLYP functionals with combinations of basis sets or using the doubly hybrid method (MUEs = 1.9-2.2 kcal/mol) give the overall best result.     

Multi-variable special polynomials using an operator ordering method

Using an operator ordering method for some commutative superposition operators, we introduce two new multi-variable special polynomials and their generating functions, and present some new operator identities and integral formulas involving the two special polynomials. Instead of calculating complicated partial differential, we use the special polynomials and their generating functions to concisely address the normalization, photocount distributions and Wigner distributions of several quantum states that can be realized physically, the results of which provide real convenience for further investigating the properties and applications of these states.     

A new non-perturbative approach in quantum mechanics for time-independent Schr?dinger equations

A new non-perturbative approach is proposed to solve time-independent Schr?dinger equations in quantum mechanics.It is based on the homotopy analysis method(HAM)that was developed by the author in 1992 for highly nonlinear equations and has been widely applied in many fields.Unlike perturbative methods,this HAM-based approach has nothing to do with small/large physical parameters.Besides,convergent series solution can be obtained even if the disturbance is far from the known status.A nonlinear harmonic oscillator is used as an example to illustrate the validity of this approach for disturbances that might be one thousand times larger than the possible superior limit of the perturbative approach.This HAM-based approach could provide us rigorous theoretical results in quantum mechanics,which can be directly compared with experimental data.Obviously,this is of great benefit not only for improving the accuracy of experimental measurements but also for validating physical theories.     

Flow physics of normal and abnormal bioprosthetic aortic valves

Flow physics of transvalvular flows in the aorta with bioprosthetic valves are investigated using computational modelling. For the efficient simulations of flow-structure-interaction in transvalvular flows, a simplified, reduced degree of freedom valve model is employed with a sharp interface immersed boundary based incompressible flow solver. Simulations are performed for normal as well as abnormal valves with reduced leaflet motion that models the effect of early leaflet thrombosis. The structure of the aortic jet and the hemodynamic stresses on the aortic wall are analysed to understand the hemodynamic impacts and possible long-term clinical implications of sub-clinical, reduced leaflet motion. The simulation results have shown that the reduced leaflet motion tilts the direction of aortic jet and generates stronger flow separation and re-attachment on the aortic wall downstream from the reduced motion leaflets. The modified flow pattern increases the wall pressure fluctuation and average wall shear stress on the downstream aortic wall, and results in the asymmetric oscillatory shear index distributions, which may have long-term clinical implications such as aortic wall damage and remodelling.     

Investigation on the effect of density ratio on the convergence behavior of partitioned method for fluid–structure interaction simulation

In order to investigate the effect of density ratio of fluid and solid on the convergence behavior of partitioned FSI algorithm, three strong-coupling partitioned algorithms (fixed-point method with a constant under-relaxation parameter, Aitken’s method and Quasi-Newton inverse least squares (QN-ILS) method) have been considered in the context of finite element method. We have employed the incompressible Navier–Stokes equations for a Newtonian fluid domain and the total Lagrangian formulation for a non-linear motion of solid domain. Linear-elastic (hyper-elastic) model has been employed for solid material with small (large) deformation. A pulsatile inlet-flow interacting with a 2D circular channel of linear-elastic material and a pressure wave propagation in a 3D flexible vessel have been simulated. Both linear-elastic and hyper-elastic (Mooney–Rivlin) models have been adopted for the 3D flexible vessel. From the present numerical experiments, we have found that QN-ILS outperforms the others leading to a robust convergence regardless of the density ratio for both linear-elastic and hyper-elastic models. On the other hand, the performances of the fixed-point method with a constant under-relaxation parameter and the Aitken’s method depend strongly on the density ratio, relaxation parameter selected for coupling iteration, and degree of deformation. Although the QN-ILS of this work is still slower than a monolithic method for serial computation, it has an advantage of easier parallelization due to the modularity of the partitioned FSI algorithm.